Experimental determination of correlations between spontaneously formed vortices in 

a superconductor 
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We have imaged spontaneously created arrays of vortices (magnetic flux quanta), generated in a 
superconducting film quenched through its transition temperature at rates around lO^K/s. From 
these images, we calculated the positional correlation functions for two vortices and for 3 vortices. 
We compared our results with simulations of the time dependent Ginzburg Landau equation in 2D. 
The results are in agreement with the Kibble-Zurek scenario of spontaneous vortex creation. In 
addition, the correlation functions are insensitive to the presence of a gauge field. 



Vortices (magnetic flux quanta) are topological defects 
of the order parameter in a superconductor. These de- 
fects are predicted to appear spontaneously, as a result 
of a rapid quench through the phase transition into the 
superconducting state. During a rapid quench, close to 
the transition temperature (Tc), the relaxation time of 
the system becomes larger then the quench time. Under 
these conditions, the system is necessarily driven out of 
equilibrium. One model that describes the outcome of 
such phase transition is the Kibble-Zurek scenario, first 
suggested by Kibble^i^ in a cosmological context. Among 
his other contributions, Zurek^ proposed terrestrial tests 
of this model in condensed matter systems, where the 
broken symmetry is U(l), such as superfluids, BEG and 
superconductors. The behavior of these systems can be 
described by a complex order parameter. Above the crit- 
ical temperature Tc, the order parameter fluctuates with 
a characteristic size of the fluctuation being ^. Is the sam- 
ple is cooled infinitely slowly towards Tc, ^ will grow until 
it reaches the size of the sample. At finite cooling rates 
the fluctuations "freeze" at some point, forming isolated, 
uncorrclated domains of the ordered state. The typical 
size of such a domain, ^ , inside which the emerging order 
parameter is coherent, depends on the cooling rate. In 
equilibrium, the order parameter in the final state should 
be uniform across the system. However, the initial mis- 
match of the phase of the order parameter between dif- 
ferent regions leads to the appearance of topological de- 
fects. In superconductors, these are vortices carrying a 
quantum of magnetic flux <I>o = h/2e. 

In the KZ model, below Tc uncorrelated domains of the 
ordered state are formed. Each of these domains picks up 
a random value of phase of the order parameter. Single- 
valuedness of the order parameter requires that the in- 
tegral of the phase accumulated along the circumference 
of any loop must be an integer multiple of 2tt. If the 
phases of different domains are random, there is a fi- 
nite probability that the phase accumulated along a loop 
around a vertex between 3 domains will be ±2tt, lead- 
ing to formation of topological defect. To calculate the 
probability of this formation, the geodesic rule is usually 



implemented. It assumes that minimal phase gradient 
will always be favorable due to minimal energy consid- 
eration. Under this assumption, only one vortex can be 
created at a vertex between 3 ordered regions. Hence, 
the geodesic rule restricts the number of topological de- 
fects created by the system. This description is limited 
to relaxation of the phase gradients and does not involve 
any additional dynamics. One of the predictions of this 
model is strong, short range correlation between vortices 
and anti- vortices. If the size of frozen fluctuations ^ is 
assumed to have a gaussian distribution, one can calcu- 
late the vortex- vortex correlation function'*. Hence, the 
fluctuations distribution above the critical temperature 
leaves its mark on the emerging vortex array. By mea- 
suring the positions of the vortices and calculating the 
correlation function it is possible to investigate the fluc- 
tuation distribution above T^. 

We note that the KZ model does not address the 
critical coarsening process which may occur after the 
transition^. In many systems, coarsening affects dras- 
tically the outcome of the transition due to defect- 
antidefect annihilation. This may not be the case for 
superconducting fllms. If the quench is fast enough, the 
system reaches low temperatures before vortices can tra- 
verse the distance to a nearby antivortex and annihilate. 
At temperatures much lower than Tc vortices are strongly 
pinned, so that any motion is practically impossible. As a 
result, vortex annihilation and coarsening is suppressed. 
Superconductors therefore offer us an unique opportunity 
to investigate the order parameter fluctuations above Tc, 
by measuring the vortex distribution after a quench. 

The validity of KZ mechanism in systems with local 
gauge symmetries (such as superconductors), has been 
questioned by several authors^. For these systems an al- 
ternative mechanism of flux trapping was suggested^Ti^. 
In this mechanism thermal fluctuations of the magnetic 
fleld are frozen inside the superconductor during the tran- 
sition. As a consequence vortices are formed in clusters of 
equal sign. The resulting vortex- vortex correlation func- 
tion should decay as a power lawiS. However, the am- 
plitude of trapped magnetic fleld fluctuations in conven- 
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tional superconducting films should be so low that vortex 
formation should be better described by the KibbleZurek 
mechanism^. 

It was suggested by several authors that topological de- 
fects in first order phase transitions^"— and in sustems 
showing spinodal decompositioni^ may form due to dy- 
namics. In this approach, the geodesic rule does not 
hold anymore. Pairs of vortices and anti-vortices are 
formed during collisions between domain walls. Although 
most of the simulations were done for first order phase 
transitions, this mechanism is claimed to be generally 
applicablei^. The correlation function for this mecha- 
nism was never calculated, but it should have two char- 
acteristic lengths: the separation between nearby vortices 
produced by domain wall collisions, and the domain size. 

Vortex pairs of opposite sign could also arise from a 
Kosterlitz-Thouless type of transitionii. In this theory, 
unbound vortex pairs appear above Tkt- If the system 
is cooled through TKT^ these vortices effectively annihi- 
late. However, if the quench is fast, some of the unbound 
vortex pairs can survive the quench, become pinned at 
low temperatures and observed. In this theory, the den- 
sity of unbound vortex pairs above Tkt increases. Con- 
sequently, the observed vortex density should depend 
on the temperature from which the system is quenched. 
Within our resolution, we found no such dependence in 
the experiment described here. 

The Kibble-Zurek (KZ) model has been tested in liq- 
uid heliumi^ii^, liquid crystals^^, in superconductors^l, 
Josephson junction a^^i^^ and superconducting loops^iS^. 
Spontaneously generated topological defects were de- 
tected in several of the aforementioned experiments. 
More sensitive testing of the model involves the deter- 
mination of correlations between these defects. There 
are only two such experiments performed to date. In 
one experiments*^, done with liquid crystals, an array of 
topological defects was imaged. However, the amount of 
data was insufficient to detect correlations beyond near- 
est neighbors. In our experiments^, we imaged sponta- 
neously formed arrays of vortices in a superconductor. 
The amount of data gathered was large, allowing a pre- 
cise determination of the two point correlations. In this 
work, we extend this study to look at correlations be- 
tween 3 defects. 

Our experimental technique is described in detail in 
previous publications^^'^^. Briefly, we image the mag- 
netic field on the surface of a superconducting film us- 
ing high resolution Magneto-optics. The superconductor 
sample consists of a 200nm thick Niobium film with Tc of 
8.9i^. The film is patterned into small squares of 
across. On top of the Nb film, we deposited a 40nm layer 
of EuSe which serves as the Magneto-Optic sensor. To 
minimize the effect of stray magnetic fields, our appara- 
tus was carefully shielded using /i-metal. The residual 
field was less than 10~^T). 

From our previous experiments^i, we know that ex- 
tremely high cooling rates are essential for spontaneous 
generation of a measurable amount of vortices. No less 




FIG. 1. Spontaneously created vortices in a superconduc- 
tor, a) A typical image of magnetic field created after a 
quench. The intensity is proportional to the local magnetic 
field. Bright and dark spots are vortices and anti-vortices 
respectively. The scale bar represents 10 /im. b) A typical 
result of the simulation of TDGL equations. The scale bar 
represents 10 ^o- 



important, fast cooling to low temperatures (far below 
Tc) traps the vortices on pinning centers, preventing an- 
nihilation of vortices and anti- vortices. Using short laser 
pulses^^, we achieved cooling rates as high as 2 • IQ'^K/s. 

Fig. la shows a typical image of spontaneously gener- 
ated vortices. The average asymmetry between the den- 
sity of positive and negative vortices is less than 1%. The 
average density of vortices was 6 ■ IO^cto"^ for a cooling 
rate of 4 • l(fK/s and 1.3 • IO^cto"^ for a cooling rate of 
2 • lO^A'/s. The scaling of the density with the cooling 
rate is consistent with the KZ model at 2D (proportional 
to the square root of the cooling rate). 

Our experimental results are compared with numer- 
ical simulations of the 2D time dependent Ginzburg- 
Landau equations. The simulations are described in de- 
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FIG. 2. The vortex- vortex correlation function G(r). Solid 
red triangles represent the correlation function calculated 
from experimental data. The statistical error bars are smaller 
than the point size. Open blue diamonds represent correla- 
tions calculated from the results of simulation. The solid line 
is a fit to theory of Liu and Mazenko^. The negative peak 
at short distance reflects vortex- ant i vortex correlations pre- 
dicted by KZ model. 



tail in2^. Briefly, we used a square 200 x 200 grid, with 
the parameters k = 1, F = 1. The initial conditions, 
A = — 0,9 — 0.7, mimic a quench from tempera- 
tures far above Tc to T = 0.7Tc. The implicit Crank- 
Nicholson scheme was employed on a staggered grid with 
step in time At = 0.01 and in space A; = 0.2. The 
time of each run was 20rGL- After this time the vor- 
tices are well defined and assumed to be pinned. A typ- 
ical result is presented in fig. lb. Bright and dark spots 
mark the positions of vortices and anti- vortices. Several 
of the parameters of the simulation are different from the 
experimental parameters. First, the cooling rate in the 
simulation is infinite, in contrast to the finite rate in the 
experiment. Finite cooling rates in the simulation scale 
the correlation length but do not affect the distribution. 
Second, the simulation is fully two dimensional, while the 
sample used in the experiment is a thin film. The inter- 
action between vortices in thin films is much stronger, 
and mediated through the magnetic field of the vortex 
outside of the film^S. 

We used images like those shown in Fig.l to deter- 
mine the correlations between the vortices. The two 
particle vortex-vortex correlation function is defined as 
G{r — r') —< n{r)n{r') > , with n{r) = 1 at the loca- 
tion of a positive vortex, —1 at the location of negative 
vortex and elsewhere. The correlation function cal- 
culated from our data is shown in Fig. 2. The distance 
in the figure was scaled by the mean vortex separation. 
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as proposed by—. Vav is related to by 



Vav = where p is the average number of vortices per 

domain. For our high cooling rate of 2 • lO^K/s, Vav is 
8.2/im. The correlation function was averaged over 260 
images, with 50,000 vortices in total. At the distance 
corresponding to the nearest neighbors the correlation 




FIG. 3. Schematic description of spontaneous formation of 
vortices on intersections between ordered domains, a) Ran- 
dom domain array, b) Close packed domain structure. Vor- 
tices and anti-vortices are marked by -|- and - respectively. 



function has a minimum. This is consistent with the 
KZ model which predicts that nearest neighbor vortices 
should have opposite polarities. The correlation function 
calculated from the simulations in the Fig. 2 was aver- 
aged over 200 realizations. Since the same scaling was 
used, the experimental and simulated correlation func- 
tions can be compared without additional parameters. 
The dashed line is a fit to the theoretical predictions by 
Liu and Mazenkc^ with ^ = 0.35ra„. In their work, Liu 
and Mazenko assumed a Gaussian distribution of the fluc- 
tuations. Even through the simulation parameters differ 
from our experimental conditions (inflnite vs. flnite cool- 
ing rate, different initial temperatures, different interac- 
tions between vortices), the resulting correlation func- 
tion is the same. Qualitatively similar correlation func- 
tion was also calculated from simulations of quenched 2D 
XY model^. The agreement with the correlation func- 
tion calculated by Liu and Mazenko"* suggests that in all 
cases, the shape of the correlation function is dictated by 
Gaussian fluctuations. 

In addition to the two point correlation function G{r — 
r'), we used our data to calculate an angular ( 3 point) 
correlation function. This function involves 3 vortices. 
For each vortex in the array, we look for its two near- 
est neighbors. Nearest neighbors are deflned as vortices 
located within a distance of ^ and have an opposite 
polarity (see fig. 3). If such neighbors are found, we cal- 
culate the angle between the two lines connecting the 
vortex to its neighbors. In the limit where the domains 
are mono-dispersed and close packed (fig. 3b), the angle 
between nearest neighbors will be 120°. For randomly 
distributed domains (fig. 3a) the angles will have a broad 
distribution. Even through the connection between the 
fluctuation distribution and the angle distribution was 
never calculated, the first obviously determines the sec- 
ond. Therefore the comparison between the angle dis- 
tribution calculated out of experimental data, and from 
the results of a simulation can be used as a validity test 
for a model. The angle distributions calculated from ex- 
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the error bars, these two distributions are consistent. 

In conclusion, we calculated the vortex-vortex correla- 
^ tion function and the angular distribution function using 

^ i " ^ ^ experimental data and compared them with the results of 

J i ^ i TDGL simulations. We found an agreement between the 

^ experimental and simulated results. This implies that 

i within experimental accuracy, the TDGL equation in 2D 

^ describes the vortex formation in our system. One conse- 

quence is that the coupling to the gauge field outside the 
sample does not affect the vortex distribution, at least 
for the parameters used in the experiment. The correla- 
tion function calculated using a quenched XY model^ 
gave qualitatively similar results, which suggests that 

gauge fields have no effect on the system evolution during 
quench. 
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FIG. 4. The angular distribution function. Solid black 
squares represent the angular distribution calculated from ex- 
perimental data. Red pentagons are angular distribution cal- 
culated from the results of the simulation. 

perimental data and from simulation are shown in Fig. 
4. Both distributions have a minimum at 9 = 0° and in- 
crease until 9 w 100° were they reach a plateau. Within 
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